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We investigate the motion of pedestrians through obscure corridors where 
the lack of visibility (due to smoke, fog, darkness, etc.) hides the precise 
position of the exits. We focus our attention on a set of basic mechanisms, 
which we assume to be governing the dynamics at the individual level. Using 
a lattice model, we explore the effects of non-exclusion on the overall exit 
flux (evacuation rate). More precisely, we study the effect of the buddying 
threshold (of no-exclusion per site) on the dynamics of the crowd and in- 
vestigate to which extent our model confirms the following pattern revealed 
by investigations on real emergencies: If the evacuees tend to cooperate and 
act altruistically, then their collective action tends to favor the occurrence of 
S ' disasters. 
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1. Introduction 

Efficient evacuation of humans from high-risk zones is a very important 
societal issue [l|. The topic is very well studied by large communities of 
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scientists ranging from logistics and transportation, civil and fire engineering, 
to theoretical physics or applied mathematics. Many modeling approaches 
(deterministic or stochastic) succeed to capture qualitatively basic behaviors 
of humans (here referred to as pedestrians ) walking within a given geometry 
towards a priori prescribed exits. For instance, social force/social velocity 
crowd dynamics models (cf. e.g. [2|, [3j, J4|), simple asymmetric exclusion 
models (see chapters 3 and 4 from JpJ as well as references cited therein), 
cellular automaton- type models p, [Zj, etco are quite efficient to evacuate 
people when exits are visible. 

The driving question of our research is: What to do when exits are not 
visible ? A search through the existing literature shows that there are studies 
done [especially for fire evacuation scenarios] on how information and way 
finding systems are perceived by individuals. One of the main questions in 
fire safety research is whether green flashing lights can influence the evacua- 



tion (particularly, the exit choice); see e.g. |8|, |9|, [l0| (and the fire engineering 
references cited therein) and [ll| (partial visibility due to a non-uniform 
smoke concentration) [12| (partial visibility as a function of smoke's temper- 
ature), [13J (flow heterogeneity due to fire spreading). On the other hand, as 
far as we are aware of, there is no literature on dynamics of pedestrians (or 
groups of pedestrians) in regions with no visibility. 

In this paper, we aim at investigating the motion of pedestrians through 
obscure corridors where the lack of visibility (due to smoke, fog, darkness, 
etc.) hides the precise position of the exits. We focus our attention on a set 
of basic mechanisms, which we assume to be governing the dynamics at in- 
dividual level. Using a lattice model, we explore the effects of non-exclusion 
on the overall exit flux (evacuation rate). More precisely, we study the effect 
of the buddying threshold (of no-exclusion per site) on the dynamics of the 
crowd and wish to investigate to which extent our model confirms the follow- 
ing pattern revealed by investigations on real emergencies: If the evacuees 
tend to cooperate and act altruistically, then their collective action tends to 
favor the occurrence of disasters^ 

By means of a minimal model, we wish to describe how a bunch of people 
located inside a dark (smoky, foggy, etc.) corridor exits through an invisible 



^Tor a detailed classification of pedestrian models, see [5j- 

2 Note that due to the lack of visibility anticipation effects (see 14|) and drifts (see [7]) 
are expected to play no role in evacuation. 



door open in one of the four walls. A few central driving questions need to 
be mentioned at this point: 

(Ql) How do pedestrians choose their path and speed when they are about to 
move through regions with no visibility? 

(Q2) Is group formation (e.g. buddying) the right strategy to move through 
such uncomfortable zones able to ensure exiting within a reasonable 
time? 

It is worth noting that this is one of the situations where group psychology 



(compare e.g. [15|, |16| and [17] ) is not really helping, especially due to the 
lack of experimental observations. On the other hand, the concept of group 
we look at here is not similar to the swarming patterns typically emphasized 
in nature by fish and or birds communities (see e.g. the 4-groups taxonomy 
in [18j, namely swarm, torus, dynamic parallel groups, and highly parallel 



groups). On top of this, note that concepts like leadership (see e.g. (19|) or 
motion fluidization by favorizing lanes formation (see e.g. [20]]) are simply 
not applicable in this context. 

The reminder of the paper has the following structure: In section [2] we 
introduce our model. Section H] presents our simulation results, while we 
close with comments, discussions and suggestions for further investigations 
in Section [5j 

2. A lattice model for the reverse mosca cieca game 

We make here first attempts in answering the questions (Ql) and (Q2). 
To do this end we employ a minimal lattice model, named the reverse mosca 
cieca gameQ, where we incorporate a few basic rules for the pedestrians mo- 
tion. 



3 Mosca cieca means in Italian blind fly. It is the Italian name of a traditional children's 
game also known as blind man's buff or blind man's bluff. The game is played in a spacious 
free of dangers area, such as a large room, in which one player, the "mosca", is blindfolded 
and moves around attempting to catch the other players without being able to see them, 
while the other players try to avoid him, hiding in plain sight and making fun of the 
"mosca" to make him/her change direction. When one of the player is finally cautgh, the 
"mosca" has to identify him by touching is face and if the person is correctly identified 
he becomes the "mosca" . Our model is a sort of reverse mosca cieca game since all the 
players (the pedestrians) group around, as if they were blindfolded, trying to catch the 
(invisible) exit. 



2.1. Basic assumptions on the pedestrians motion 

We assume that the motion of pedestrians is governed by the following four 
mechanisms: 

(Al) In the core of the corridor, people move freely without constraints; 

(A2) The boundary is reflecting, possibly attracting; 

(A3) People are attracted by bunches of other people up to a threshold; 

(A4) People are blind in the sense that there is no drift (desired velocity) 
leading them towards the exit. 

(Al)-(A4) intend to describe the following situation: 

Since, in this framework, neighbors (both individuals or groups) can not 
be visually identified by the individuals in motion, basic mechanisms like 
attraction to a group, tendency to align, or social repulsion are negligible 
and individuals have to live with "preferences" . Essentially, they move freely 
inside the corridor but they like to buddy to people they accidentally meet 
at a certain point (site). The more people are localized at a certain site, the 
stronger the preference to attach to it. However if the number of people at a 
site reaches a threshold, then such site becomes not attracting for eventually 
new incomers. We refer to (A3) as the buddying mechanism. 

Once an individual touches a wall, he/she simply felts the need to stick 
to it at least for a while, i.e. until he/she can attach to an interesting site 
(having conveniently many hosts) or to a group of unevenly occupied sites 
or the exact location of the door is detected. 

Since people have no desired velocity, their diffusion (random walk) to- 
gether with the buddying are the only transport mechanisms. Can these 
eventually lead to evacuation? 

In the following, we study the effect of the threshold (of no-exclusion 
per site) on the overall dynamics of the crowd. We will describe our results 
in terms of macroscopic quantities like averaged outgoing flux, stationary 
occupation numbers, and stationary correlations. 

2.2. The lattice model 

To describe our model, we start off with the construction of the lattice. It is 
worth mentioning that the approach we take here is very much influenced by 
a basic scenario described in [21J for randomly moving sodium ions willing 



to pass through a switching on-off membrane gate. The major difference 
here is twofold: the gate is permanently open and the buddying principle is 
activated. 

Let ei := (1,0) and €2 = (0, 1) denote the coordinate vectors in M 2 . Let 
A C Z 2 be a finite square with odd side length L. We refer to this as the 
corridor. Each element a; of A will be called a cell or site. The external 
boundary of the corridor is made of four segments made of L cells each; the 
point at the center of one of these four sides is called exit. 

Let N be positive integer denoting the (total) number of individuals inside 
the corridor A. We consider the state space X := {0, . . . , N} A . For any state 
n G X, we let n(x) be the number of individuals at cell a|j. 

We define a Markov chain n t on the finite state space X with discrete 
time t = 0, 1, . . .. The parameters of the process will be the integers (possibly 
equal to zero) Q,T, W > called respectively minimal quantum, threshold, 
and wall interaction, and the real number R G [0, 1] called the rest parameter. 
We finally define the function S : N — > N such that 

s(k)-l Q [ik>T 

for any fceN. Note that for k = we have S(0) = Q. 

At each time t, the A^ individuals move simultaneously within the corridor 
according to the following rule: 

For any cell x situated in the interior of the corridor A, and all y nearest 
neighbor of x, with n G X, we define the weights 

w(x,x) := RS(n(x)) and w(x,y) := S(n(y)). 

Also, we obtain the associated probabilities 

p(x,x) and p(x,y) 

by dividing the weight by the normalization 

2 2 

w(x, x) + y^w(x, x + ei) + 2_\ w{x, x — e^. 

i=l i=l 



4 Note that, at each cell x, we don't control how many individuals are located. We only 
know that n(x) <E [0,iV] for all x 6 A. 



Let now x be in one of the four corners of the corridor A, and take y as one 
of the two nearest neighbors of x inside A. For n G X, we define the weights 

w(x, x) := R(S(n(x)) + 2W) and w(x, y) := S(n(y)) + W 

and the associated probabilities p(x, x) and p(x, y) obtained by dividing the 
weight by the suitable^ normalization. 

It is worth stressing here that T is not a threshold in n(x) - the number 
of individuals per cell. It is a threshold in the probability that such a cell is 
likely to be occupied or not. 

For x G A neighboring the boundary (but neither in the corners, nor 
neighboring the exit), y one of the two nearest neighbor of x inside A and 
neighboring the boundary, z the nearest neighbor of x in the interior of A, 
and n G X, we define the weights 



w(x, x) 
w(x,y) 
w(x, z) 



R{S{n(x)) + W) 
S{n{y)) + W 
S(n(y)). 



The associated probabilities p(x, x), p(x,y), and p(x,z) are obtained by di- 
viding the weight by the suitable normalization. 

Finally, if x is the cell in A neighboring the exit and y is one of the two 
nearest neighbor of x inside A and neighboring the boundary, z being the 
nearest neighbor of x in the interior of A, and n G X, we define the weights 

w(x, x) := RS{n(x)) 
w(x, y) := S(n(y)) + W 
w(x,z) := S(n(y)) 
w(x, exit) := T + Q. 

The associated probabilities p(x, x), p(x,y), p(x,z), and p(x,exit) are ob- 
tained by dividing the weight by the suitable normalization. 

The dynamics is then defined as follows: 

At each time t, the position of all the individuals on each cell is updated 
according to the probabilities defined above. If one of the individuals jumps 
on the exit cell a new individual is put on a cell of A chosen randomly with 
the uniform probability 1/L 2 . 



5 To get a probability, we divide the weights by the sum of the involved weights: It is 
about 3 terms if the site lies in a corner of the corridor, 4 terms if it is close to a border, 
and then 5 terms for a pedestrian located in the core. 

6 



3. Interpretation of the buddying threshold T 

The possible choices for the parameter T correspond to two different 
physical situations. The first one, for T = 0, the function S(k) is equal to q 
(the minimal quantum) whatever the occupation numbers^ are. This means 
that each individual has the same probability to jump to one of its nearest 
neighbors or to stay on his site. This is resembling the independent symmetric 
random walk case; the only difference is that with the same probability the 
individuals can decide not to move. We expect that this "rest probability" 
just changes a little bit the time scale. 

The second physical case is T > 0. For instance, T = 1 means mild 
buddying, while T = 100 would express an extreme buddying. No simple 
exclusion is included in this model: on each site one can cluster as many 
particles (pedestrians) as one wants. The basic role of the threshold is the 
following: The weight associated to the jump towards the site x increases 
from Q to Q + T proportionally to the occupation number n(x) until n(x) = 
T, after that level it drops back to Q. Note that this rule is given on weights 
and not to probabilities. Therefore, if one has T particles at y and T at 
each of its nearest neighbors, then at the very end one will have that the 
probability to stay or to jump to any of the nearest neighbors is the same. 
Differences in probability are seen only if one of the five (sitting in the core) 
sites involved in the jump (or some of them) has an occupation number large 
(but smaller than the threshold). 

4. Results 

We have studied numerically the model described in section [2] in the case 
of purely reflecting boundary and eventually resting individuals, that is we 
have considered the case 

R = 1 and W = 0. 

We have also fixed to one the value of the minimal quantum and to 101 the 
length of the side of the square A, that is 

Q = 1 and L = 101. 
6 The definition of occupation numbers is given in section f4. 2 1 
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Figure 1: Averaged outgoing flux vs. time in the case T — and N — 100 on the left and 
T = 100 and N = 100 on the right. The inset is a zoom in the time interval [4 x 10 6 , 5 x 10 6 ] 
on the left and [1.4 x 10 7 , 1.5 x 10 7 ] on the right. 



We have studied the behavior of the system by varying the treshold param- 
eter, in particular we have considered the cases 

T = 0,1, 5, 30, 100. 

It is worth noting that in the case T = no buddying effect is introduced into 
the model, that is the individuals behave as standard independent random 
walkers. The number of individuals has been also varied; for each value of 
the treshold we have studied the cases 

N = 100, 600, 1000, 6000, 10000. 

For the peculiar values T = 30, 100 more cases have been taken into account. 
More precisely, for the choices T = 30 and T = 100 we have also analyzed 
the supplementary cases 



and 



respectively. 



N = 2000, 2200, 2400, 2600, 2800, 3000, 3300 
N = 1300,1600,2000,3000, 



4-1. Average outgoing flux 

The first interesting quantity that one has to compute is obviously the average 
outgoing flux that is to say the ratio between the number of individuals which 
exited the corridor in the time interval [0, £f] and if . 

This quantity fluctuates in time, but for times large enough it approaches 
a constant value. In order to observe relative fluctuations smaller than 10 -2 
we had to use tj = 5x 10 6 . To capture the extreme buddying case T = 100, 
we used if = 1.5 x 10 7 (see figured]). 
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Figure 2: Averaged outgoing flux vs. number of pedestrians. The symbols o, x, *, □, and 
+ refer respectively to the cases T — 0, 1, 5, 30, 100. The straight line has slope 8 x 10 -6 
and has been obtained by fitting the Monte Carlo data corresponding to the case T = 0. 



Figure |2] depicts our results, where the averaged outgoing flux is given 
as a function of the number of individuals. At T = 0, that is when no 
buddying between the individuals is put into the model, the outgoing flux 
results proportional to the number of pedestrians in the corridor; indeed the 
data represented by the symbol o in Figure [2] have been perfectly fitted by a 
straight line. 

The appearance of the straight line was expected in the case T = since 
in this case the dynamics reduces to that of a simple symmetric random 
walk with reflecting boundary conditions. This effect was studied rigorously 
in the one-dimensional case and via Monte Carlo simulations in dimension 



two in |21[. The order of magnitude of the slope can be guessed with a simple 
argument: in two dimension the typical time for a walker started to any point 
in A to reach the boundary is L 2 and the number of times the walker has 
to touch the boundary to reach the exit is L. So that the typical time for a 
walker to reach the exit is L 3 . This argument yields that the outgoing flux is 
of order N/L 3 . Since in the simulation we used L = 101, we get an estimate 
of the slope approximatively eight time smaller with respect to the measured 
one. 

When a weak buddying effect is introduced in the model, that is in the 
case T = 1, we find that if the number of individuals is small enough, say iV < 
6000, the behavior is similar to the one measured in the absence of buddying 
(T = 0). At N = 10000, on the other hand, we measure a larger flux; 
meaning that in the crowded regime small buddying favors the evacuation of 
the corridor [i.e. it favors the finding of the door]. 

The picture changes completely when buddying is increased. To this end, 
see the cases T = 5, 30, 100. The outgoing flux is slightly favored when the 
number of individuals is low and strongly depressed when it this becomes 
high. The value of iV at which this behavior changes strongly depends on 
the threshold parameter T. 

4-2. Stationary occupation number 

In order to have a deeper insight in the behavior of the system we have 
computed the so called stationary occupation numbers. We let (•) be the 
stationary average of a random variable associated with the system and u(x), 
for any x G A, the random variable giving the number of particles n(x) on 
the site x divided times the density N/L 2 (with this normalization we expect 
to get typical values close to one). We define stationary occupation number 
at site x as the stationary mean (u(x)). 

From the computational point of view this task is not an easy one for two 
different reasons: 

The presence of the exit in the middle of one of the four boundaries 
makes the system not translationally invariant. Therefore, we expect 
that the occupation number (u(x)) does depend on the site x; 

- Computing {u(x)) amounts to give a Monte Carlo estimate of the sta- 
tionary measure of the Markov Chain. 



10 



In order to estimate the occupation number, we have run jobs of length 
if = 3 x 10 7 and measured the occupation number each tmeas — 100 Monte 
Carlo steps after having waited tterm = 100000 to let the system termalize 
and loose the memory of the starting condition. 

The parameter tmeas was chosen equal to 100 after analyzing the au- 
tocorrelation functions. Given a sequence of data m(j) with j = 1, . . . , J 
obtained at times j, the autocorrelation function [22] is defined as 



j , j 



for £ = 0, . . . , J' with J' <C J, where Var is the variance of the sample defined 

as 

j j 

Var := - ^{m(i)f - \j^m{j) 

i=i i=i 

Autocorrelations are supposed to decay exponentially, so that the autocor- 
relation time is defined as the smallest time such that the autocorrelation is 
smaller than 1/e. When stationary observables are computed via a Monte 
Carlo scheme, one constructs averages by sampling data separated in time 
by at least one correlation time, preferably two. 

In figure[3]we have depicted the autocorrelation of the occupation number 
computed at the center of the lattice, at the four sites on the axes parallel 
to the coordinates ones and passing through the center of the corridor at 
distance [L/4\ from the center itself, and at the four sites on the axes parallel 
to the coordinates ones and passing through the center of the corridor at 
distance |_£/2J from the center itself. We have depicted the data in the cases 
T = 5, 30, 100 and iV = 1000, 10000. In all the cases the correlation time 
is smaller than 100. So that chosing tmeas = 100 seems to be a reasonable 
choice. 

We have computed the stationary occupation number in the center of the 
corridor A, that is at the site whose coordinates are both equal to |_-^/2j + 
1, where \a\ denotes the largest integer smaller than the real number a. 
Recall that L is an odd number. Morover we have computed the stationary 
occupation number associated with all the sites of the corridor along the two 
straight lines parallel to the coordinate axes and passing through the center 
of the corridor itself. 



11 
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Figure 3: Autocorrelation vs. time for the occupation number at the sites described in the 
text. Plots refer to the case N = 1000, 10000 from the left to the right and T = 5, 30, 100 
from the bottom to the top. Solid, dashed, and dotted lines refer respectively to the center 
of the lattice, to sites at distance \_L/4\ from the center, and to sites at distance [L/2\ 
from the center. 
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Figure 4: Stationary occupation number vs. distance from the center. The symbols +, 
x, o, and □ refer respectively to upward, downward, leftward, and rightward sites with 
respect to the center of the lattice. Plots refer to the case N = 100, 1000, 10000 from the 
left to the right and T — 0, 1, 5 from the bottom to the top. 
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Figure 5: Stationary occupation number vs. distance from the center. The symbols +, 
x, o, and □ refer respectively to upward, downward, leftward, and rightward sites with 
respect to the center of the lattice. Plots refer to the case N = 100, 1000, 10000 from the 
left to the right and T = 30, 100 from the bottom to the top. 



14 



Data in Figures @] and Figure \5\ are plotted as a function of the distance 
from the center of the corridor. We have reported on the horizontal and 
vertical axes the distance of the site from the center of the corridor and the 
corresponding stationary occupation number, respectively. The symbols +, 
x, o, and □ refer respectively to upward, downward, leftward, and rightward 
sites with respect to the center of the lattice. The exit of the corridor is to 
the left of the center of the corridor at distance \_L/2\ + 1 = 51. 

In all the pictures, we note that the upward and downward behaviors are 
similar; indeed we do not expect any vertical asymmetry. On the other and 
the rightward and leftward behaviors are different with respect to each other 
and, also, with respect to the vertical behavior. This is quite intuitive since 
the exit has been placed on the horizontal line passing through the center of 
the corridor. 

In all the pictures it is clearly present a depletion effect in the rightward 
direction. This was obviously expected. This effect does not depend on N, 
while it strongly depends on T. Particularly, we note that this gets more 
and more important with increasing T. It is also worth nothing that the 
horizontal leftward behavior departs from the vertical one at large threshold 
T. 

Fixed T = 0, 1, 5, the three corresponding graphs at iV = 100, 1000, 10000 
are very similar among each others, see figure HJ On the other hand at 
T = 30, 100 the behavior at N = 10000 is somehow singular, see figure |5J 
This suggests that for large threshold T and number of individuals N the 
system reaches a stationary state different from the one characterizing the 
small threshold and/or small number of pedestrian regimes. This remark 
sounds reasonable also in view of the outgoing flux behaviors recorded in 
figure [2J 

A possible explanation for this regime is that pile of individuals are formed 
and the outgoing flux is controlled by the pile dynamics, which is typically 
much slower than the single partcle one. 

4-3. Stationary correlations 

Another relevant quantity we have studied at equilibrium is the correlation 
between the occupation number random variables. Recall that (•) denotes 
the stationary average of a random variable associated with the system. For 
any x, y G A we call truncated correlation of the occupation numbers at sites 
x and y the quantity 

([u(x) - (u(x))][u(y) - (u(y))}} = (u(x)u(y)} - (u(x))(u(y)} 
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Figure 6: Stationary truncated correlations vs. distance from the center. The symbols +, 
x, o, and □ refer respectively to upward, downward, leftward, and rightward sites with 
respect to the center of the lattice. Plots refer to the case N = 100, 1000, 10000 from the 
left to the right and T = 0,5, 100 from the bottom to the top. 
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that is to say the covariance between those random variables. Note that for 
x = y we simply get the variance of the random variable u(x). 

In Figure |6j we have plotted the normalized truncated correlations of 
the occupation numbers of the site in the center of the corridor and a site 
on the vertical and on the horizontal straight line parallel to the coordinate 
axes passing through the center of the corridor itself. In all the figures, 
the truncated correlations have been normalized times the covariance of the 
occupation number at the center of the corridor. 

The information in Figure |6] is represented as a function of the distance 
from the center of the corridor. We show there on the horizontal and vertical 
axes the distance of the site from the center of the corridor and the corre- 
sponding normalized truncated correlation, respectively. The symbols +, x, 
o, and □ refer respectively to upward, downward, leftward, and rightward 
sites with respect to the center of the lattice. 

The relevant remark at this point is that at T = no correlation is 
detected, while at T > nonzero correlations start to occur at sufficiently 
large N. 

5. Discussions 

In this section, we only wish to emphasize a few basic aspects concerning 
the dynamics of pedestrians within regions with no visibility that we discover 
with our modeling of buddying: 



- In [23], the authors present an experiment whose purpose was to study 
evacuees exit selection under different behavioral objectives. The ex- 
periment was conducted in a corridor with two exits located asym- 
metrically. This geometry was used to make most participants face a 
nontrivial decision on which exit to use. The statistical approach used 



m 



231 ] seems to indicate that the individuals of an evacuating crowd 
may not be able to make optimal decisions when assessing the fastest 
exit to evacuate. In their case, the evacuation (egress) time of the 
whole crowd turned out to be shorter when the evacuees behave egois- 
tically instead of behaving cooperatively. This is rather intriguing and 
counter intuitive fact, and it is very much in the spirit of the effect of 
the threshold T we observed in section HJ 

- Note that for low densities the buddying mechanism increase the out- 
going flux: compare points and straight line in the pictures (straight 
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line is essentially the not-buddying case), while at large densities the 
scenario is dramatic: isolated individuals may turn to have a bigger es- 
cape chance than a large group around a leader [behavior recommended 
by standard manuals on evacuation strategies, see e.g. [24j, p. 122.]. 
This suggests that evacuation strategies should not rely too much on 
the presence of a leader. 

- Quite probably, in order to decide what is the best strategy for a given 
scenario cooperation (grouping, buddying, etc.) or selfishness (walking 
away from groups) one would need to explore also basic aspects con- 
nected to the prisoner's dilemma, or to a broader context like stochastic 
game theory |25| . 
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